function z = hwozd_zmi(I,n,m)
    addpath('..\kombFactSnms');
    m = abs(m);
    z = 0;
    k = (n - m)/2;
    for s = 0:1:k
        b = (n - m)/2 - s;
        for a = 0:1:b
            for d = 0:1:m
                num1 = komb(m,d) * komb(b,a) * ((-1)^d);
                num2 = hwozd_Snms(n,m,s) * hwozd_cm(I, n - 2*s - 2*a - d, 2*a + d);
                z = z +  num1 * num2;
            end;
        end;
    end;
    z = z * (n + 1)/ pi;
